'''
This class represents grand canonical ensemble representation of the system
'''
import hamiltonian

class GrandCanonical:
    def __init__(self, kMatrix, u, mu):
        self.nSites = kMatrix.nSites
        self.measurements = []
        for nElectrons in range(2 * self.nSites + 1):
            h = hamiltonian.Hamiltonian(kMatrix, uList, mu, nElectrons)
            self.measurements += [h.meauserments]